Probing Phase Fluctuations in a 2D Degenerate Bose Gas by Free Expansion 
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We measure the power spectrum of the density distribution of a freely expanding 2D degenerate 
Bose gas, where irregular density modulations gradually develop due to the initial phase fluctuations 
in the sample. The spectrum has an oscillatory shape, where the peak positions are found to be 
independent of temperature and show scaling behavior in the course of expansion. The relative 
intensity of phase fluctuations is estimated from the normalized spectral peak strength and observed 
to decrease at lower temperatures, confirming the thermal nature of the phase fluctuations. We 
investigate the relaxation dynamics of nonequilibrium states using the power spectrum. Free vortices 
are observed with ring-shaped density ripples in a perturbed sample after a long relaxation time. 
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Phase coherence is one of the main characteristics of 
superfluidity and critically affected by the dimensions of 
the system. In low dimensional systems, large thermal 
and quantum phase fluctuations prohibit the establish- 
ment of long-range phase coherence which is a typical 
order parameter for a three-dimensional superfluid [TJ [2] . 
Nevertheless, a two-dimensional (2D) interacting system 
can undergo a superfluid phase transition at a finite crit- 
ical temperature with an algebraically decaying coher- 
ence. This transition is successfully described in the 
Berezinskii-Kosterlitz-Thouless (BKT) theory [3, 4 as a 
topological phase transition, where the critical point is 
associated with spontaneous pairing of free vortices with 
opposite circulations. 

Recent experiments with 2D atomic Bose gases have 
demonstrated that this BKT physics can be studied in 
a finite-size trapped sample [5HTT] . Phase coherence and 
thermodynamic properties have been investigated using 
matter-wave interference and by detailed analysis of the 
in situ density and momentum distributions of trapped 
samples, observing the algebraic decay of coherence [5], 
a presuperfluid regime [7HT0] , and the scale invariance of 
the equation of state [9] [11] . It is now highly desirable 
to have quantitative probes directly sensitive to phase 
fluctuations in order to study the topological nature of 
the phase transition. In particular, nonequilibrium phase 
dynamics near the critical point would provide valuable 
insights on the BKT transition. Relaxation dynamics in 
the 2D XY model have been under intense theoretical in- 
vestigation p~2j [13] and recently an experimental scheme 
to study a dynamic BKT transition in 2D Bose gases was 
proposed [T4] . 

In this paper, we demonstrate a new quantitative probe 
for phase fluctuations in a 2D degenerate Bose gas using 
the density correlations in a freely expanding sample. Ir- 
regular density modulations gradually develop in the co- 
herent part of the sample during expansion. We observe 
that the power spectrum of the density distribution has 
an oscillatory shape and find that the peak positions are 



independent of temperature and show scaling behavior in 
the course of expansion. The relative intensity of phase 
fluctuations is estimated from the spectral peak strength 
normalized with the central density of the coherent part 
and we show that it decreases at lower temperature in 
thermal equilibrium. This confirms the thermal nature 
of phase fluctuations in the 2D system. In addition, 
we investigate the relaxation dynamics of nonequilibrium 
states by measuring the time evolution of the relative in- 
tensity of phase fluctuations. Interestingly, ring-shaped 
density ripples stochastically appear in a perturbed sam- 
ple after a long relaxation time, which we identify with 
free vortices having a long lifetime in a 2D sample. 

We prepare a 2D degenerate Bose gas of 23 Na atoms in 
a single pancake-shape optical dipole trap [I5j [16] . Ther- 
mal atoms in the \F = l,mf = —1) state are loaded 
from a plugged magnetic trap [17 into the optical trap 
and evaporative cooling is applied by reducing the trap 
depth. The sample temperature is controlled by the final 
trap depth in the evaporation, resulting in 0.6 to 1.3 x 10 6 
atoms in a sample. Finally, the optical trap depth ramps 
up and the trapping frequencies (uo x1 uj yi uj z ) = 27rx(3.0, 
3.9, 370) Hz. The cooling procedure is intentionally set 
to be slow over 15 s, ensuring thermal equilibrium. The 
lifetime in the optical trap is over 50 s. In the Thomas- 
Fermi approximation, the chemical potential is about 
h x 260 Hz less than the confining energy huj z , so we ex- 
pect 2D physics in the phase coherence of the sample at 
low temperature. The dimensionless interaction strength 
g = ay/S7rmuj z /h ~ 0.013, where a is the 3D scattering 
length and m is the atomic mass. The in-plane density 
distribution n(x,y) is measured by taking an absorption 
image after an expansion time t e which is initiated by 
suddenly turning off the trapping potential. 

Expansion has been a conventional and powerful 
method in quantum gas experiments to study coher- 
ence properties of a sample [7J [10]. In our experiment, 
we are interested in the short expansion regime where 
t e <C l/uj Xy y and the phase coherence information would 
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FIG. 1: (color online) Emergence of density fluctuations in 
a freely expanding 2D Bose gas. Density distributions af- 
ter (a) t e — 14 ms and (d) 23 ms of time-of- flight. Density 
fluctuations gradually develop during expansion, increasing 
their length scale and visibility. (b,e) The horizontal density 
profiles in the center of the samples. The red dashed lines 
indicate the averaged profiles over 10 individual realizations 
of the same experiment. The coherent fraction 77 w 30% (see 
text for details). The power spectrum of the density distri- 
bution is measured with the magnitude square of its Fourier 
transform. (c,f) The averaged power spectra corresponding 
to (a,d). 



be revealed as density correlations. Since the sample ini- 
tially expands fast along the tight direction, the atom in- 
teraction effects are rapidly reduced and the subsequent 
evolution in other directions can be described as free ex- 
pansion. When the sample contains phase fluctuations, 
self interference would result in density modulations in 
the expanding sample. This method has been exploited 
in previous studies of phase fluctuations in elongated 
Bose-Einstein condensates [18] and ID Bose gases [19]. 

We observe that density fluctuations develop in an 
expanding 2D Bose gas (Fig. [I]), where the character- 
istic size and the visibility of the density lumps increases 
with the expansion time. Density fluctuations appear 
discernible only when the sample shows a bimodal den- 
sity distribution so we refer the center part as the co- 
herent part of the sample in the following [6j [7]. In 
order to obtain the density correlation information, we 
measure the power spectrum of the density distribution 
as the square of the magnitude of its Fourier transform, 
P(o) = I J dxdye % ^' r n(f)\ 2 . Although the spatial pattern 
of the density modulations appears random in each re- 
alization, the power spectrum clearly reveals a multiple 
ring structure that scales down with the expansion time 
[Fig. [TJc) and (f)]. For quantitative analysis, we obtain 
an ID spectrum P(q) by azimuthally averaging the 2D 
spectrum (Fig. [2| [20 . The strong signal around q = 
corresponds to the finite size of the coherent part. 

The scaling behavior of the spectrum can be qualita- 
tively understood in terms of the Talbot effect [2TJ [22] . It 
is well known in near-field diffraction that when a grating 
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FIG. 2: (color online) The power spectrum of density fluctua- 
tions, (a) ID power spectra P(q) are obtained by azimuthally 
averaging the averaged 2D spectra for various expansion times 
t e . Each spectrum is displayed with an offset for clarity. The 
dashed lines are the corresponding spectra for the averaged 
density distributions over 10 individual realizations, (b) P(q) 
at t e = 17 ms for various temperatures. The spectral peak 
positions g n 's are determined from a fit of multiple gaussian 
curves to each spectrum (n is the peak order number). The 
vertical dashed lines indicate the average peak positions, (c) 
The empirical function hqnt e /2nm — an 1 are fit to the peak 
positions in (a) for each t e . (d) The expansion time depen- 
dence of the two free parameters a and 7. 



is illuminated by monochromatic waves, the identical self- 
image of the grating is formed at a distance L T = 2cP/\ 
away from the grating, where d is the grating period 
and A is the wavelength of the incident radiation. The 
same effect occurs with a phase grating [23]. In matter 
wave optics, A = h/mv, where v is the incident speed 
of atoms [24] , so the propagation time for self-imaging is 
defined as t^ = L T /v = 2md 2 /h independent of v. If 
we consider a 2D Bose gas as a macroscopic matter wave 
containing phase fluctuations at all length scales, it is ex- 
pected that the component of wavenumber q satisfying 
the Talbot condition q 2 = 4:irm/ht e will emerge predom- 
inantly in the density distribution at a given expansion 
time t e . The multiple peaks in P(q) can be accounted 
for by the fractional Talbot effect where self-images with 
smaller periods d/n (n > is an integer) are produced 
at L T /2n [22]. 

Recently, theoretical calculations on the spectrum of 
density modulations have been performed for a homoge- 
neous 2D Bose gas at low temperatures [25] , showing that 
the nth peak position q n closely satisfies hq 2 l t e /27rm « 
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FIG. 3: (color online) Temperature dependence of phase 
fluctuations at thermal equilibrium. The relative intensity of 
phase fluctuations is estimated from the normalized strength 
of the first spectral peak P(qi) and the relative temperature 
is parameterized with the coherent fraction. The upper insets 
show the density profiles for different temperatures. The red 
dashed lines are gaussian fits to the thermal wings, from which 
temperatures are estimated (lower inset). Each data point 
consists of ten independent measurements and error bars in- 
dicate standard deviation. 



(n — 1/2). In particular, they predict that for sufficiently 
long expansion times the spectrum remains self-similar 
during expansion and its shape is determined only by 
the exponent of the power-law decay of the first order 
coherence function. 

In our experiment, we observe that the spectrum pre- 
serves its oscillatory shape during expansion and that 
the peak positions q n 's are independent of temperature 
(Fig.|2|, which are in qualitative agreement with the the- 
oretical prediction. However, we find different scaling 
behavior of g n 's in the measured spectra, which is well 
described as hq^t e /27rm = an 1 with 0.2 < a < 0.45 and 
0.7 < 7 < 1 for t e = 10 - 25 ms [Fig. [2jd)]. Further- 
more, the phase of the spectral oscillation is opposite 
to the theoretical prediction, suggesting that an addi- 
tional peak is present at qo = 0, which is also hinted 
by the shoulder-like hump in q < qi/y/2. We rule out 
the finite size effects by seeing no dependence of q n 's on 
sample size as well as temperature, implying that the ob- 
served scaling behavior might be intrinsic to the expan- 
sion dynamics. We note that the long expansion time 
condition ^ht e /m ^> for the validity of the theoret- 
ical prediction is marginally fulfilled in our experiment 
(^/ht e /m/^2D where = fr/y/mji is the 2D heal- 
ing length, fi being the chemical potential [26] . 

The universality of the spectral peak positions suggests 
that the spectral peak strength can be used as a measure 
of the magnitude of phase fluctuations in a sample. In 
order to quantify the relative intensity of phase fluctua- 



tions, we normalize the strength of the first spectral peak 
with the square of the central density n c of the coherent 
part in the sample, P(qi) = P{q\)/n 2 c , where n c is de- 
termined from a fit of two gaussian curves to the density 
distribution [27H29] . 

Using P(<2i), we first investigate the temperature de- 
pendence of phase fluctuations in thermal equilibrium. 
To estimate the relative temperature to the critical point 
in a model-independent way [8j [11], we use the coher- 
ent fraction 77 that is defined as the ratio of the atom 
number of the coherent part to the total atom number. 
The value of 77 was constant within 5% for our expansion 
times. Fig. [3] shows that P{q\) is suppressed at lower 
temperature (higher 77), confirming the thermal nature 
of phase fluctuations. 

For a weakly interacting 2D Bose gas, especially for our 
small g = 0.013, the BKT critical temperature T c is close 
to the BEC critical temperature for a trapped ideal Bose 
gas T c ,bec = 0Mh(Lj x LjyLj z N)^ 3 /k B [3E3GH]. We es- 
timate T c « 80 nK for N « 1.3xl0 6 . Since k B T c « 4faj z , 
thermal populations in the tight direction is not negligi- 
ble, accounting for the gaussian-like profile of the satu- 
rated thermal cloud [27l [28] . The rapid increase of P(qi) 
at 77 < 0.2 might indicate the behavior in the proxim- 
ity of the critical point. In Ref. [7 , the critical point 
was identified at 77 « 0.1 with the abrupt change in the 
width of the coherent part. Since the spatial extent of 
the coherent part becomes small, the signal-to-noise ratio 
is poor when 77 < 0.05 so we cannot study the presuper- 
fluid regime where the decay of the coherence function 
changes from algebraic to exponential, which might be 
reflected in the spectral shape. 

The power spectrum can be used to study nonequi- 
librium dynamics in a 2D Bose gas. For this study, we 
prepare a 2D sample in a nonequilibrium state by trans- 
ferring a condensate instead of thermal atoms from the 
plugged magnetic trap into the optical trap. The induced 
perturbations are small enough that the density profile in 
the optical trap is quite close to that at equilibrium. In 
Fig. [4ja), we plot the time evolutions of samples in vari- 
ous initial conditions in the plane of P(qi) and 77 (Fig.[Z]f), 
clearly showing that the nonequilibrium states decay to 
equilibrium. The decay time of the excess phase fluctua- 
tions with respect to the equilibrium value is measured to 
~ 4 s, corresponding to ~ 10 collision times in our typical 
condition. Note that the hottest sample first decays and 
then moves along the equilibrium line with increasing 77 
because of the evaporation cooling due to the finite trap 
depth. This verifies that our previous measurements are 
indeed for phase fluctuations in thermal equilibrium. 

Remarkably, ring-shaped density ripples are observed 
in the perturbed samples after long relaxation times 
(Fig. |4]). We believe that this corresponds to vortex ex- 
citations generated in the sample transferring procedure. 
Since a vortex in 2D can decay only via pairing with 
another vortex with opposite circulation or drifting out 
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FIG. 4: (color online) Relaxation of a nonequilibrium 2D 
Bose gas. (a) Time evolutions of various nonequilibrium 
states in the plane of P(qi) and rj (t e = 17 ms). The dashed 
line is the interpolation for the thermal equilibrium data with 
t e = 17 ms in Fig. [3] (b) The decay curves of the excess phase 
fluctuations P(qi) — -P eq ((7i), where P e q(qi) is the equilibrium 
value corresponding to 77 at th = 18 s of hold time. Density 
distributions of the lowest temperature sample at (c) th = s 
and (d) 18 s with t e = 17 ms, and at th = 18 s with (e) 
t e = 8 ms and (f) 23 ms. The samples in (e) and (f) were 
strongly perturbed to see more ring-shaped density ripples. 



of the finite sample, we may expect a metastable state 
with vortices having a long lifetime. Recently, it has 
been reported that vortex excitations survive longer in 
an oblate condensate because Kelvin mode excitations 
on a vortex line are suppressed [30l |3T] . The ring pattern 
appears more often at lower temperature with stronger 
perturbations. In the two lowest temperature cases (cir- 
cle and triangle in Fig. [4|, the appearance probability is 
about 60% at th = 18 s of hold time, where the spectral 
strength of the samples with vortices is about 10% higher 
than without them. In Fig. |4jb), the decay rates of the 
two coldest samples become slightly slower after th > 5 s, 
which might be attributed to the long lifetime of vortices. 
Fig. gd)-(f) shows the expansion dynamics of the ring- 
shape density ripples. In the thermal equilibrium case, 
the ring pattern was never seen for 77 > 0.2. 

In conclusion, we have demonstrated the power spec- 
trum of the density fluctuations in a freely expanding 
2D Bose gas as a new quantitative probe for phase 
fluctuations. Together with more controlled perturba- 
tions [32j [33] , we expect this method to be extended 
for studying nonequilibrium phenomena in BKT physics 
such as critical exponents [12j [13[ and dynamic transi- 
tions P3J [34] . 

This work was supported by the NRF grants funded 
by the Korea government (MEST) (Nos. 2010-0010172, 
2011-0017527, 2008-0062257, and WCU-R32- 10045). JC, 



SSW, and YS acknowledge support from the Global PhD 
Fellowship, the Kwanjeong Scholarship, and the TJ Park 
Science Fellowship, respectively. 



Electronic address: yishin@snu.ac.kr 
[1] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 
(1966). 

[2] P. Hohenberg, Phys. Rev. 158, 383 (1967). 

[3] V. Berezinskii, Sov. Phys. JETP 34, 610 (1972). 

[4] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 

(1973); J. M. Kosterlitz, ibid. 7, 1046 (1974). 
[5] Z. Hadzibabic, P. Kriiger, M. Cheneau, B. Battelier, and 

J. Dalibard, Nature (London) 441, 1118 (2006). 
[6] P. Kriiger, Z. Hadzibabic, and J. Dalibard, 

Phys. Rev. Lett. 99, 040402 (2007). 
[7] P. Clade, C. Ryu, A. Ramanathan, K. Helmerson, and 

W. D. Phillips, Phys. Rev. Lett. 102, 170401 (2009). 
[8] S. Tung, G. Lamporesi, D. Lobser, L. Xia, and E. A. Cor- 
nell, Phys. Rev. Lett. 105, 230408 (2010). 
[9] C.-L. Hung, X. Zhang, N. Gemelke, and C. Chin, Nature 
(London) 470, 236 (2011). 
[10] T. Plisson et a/., Phys. Rev. A 84, 061606(R) (2011). 
[11] T. Yefsah, R. Desbuquois, L. Chomaz, K. J. Giinter, and 

J. Dalibard, Phys. Rev. Lett. 107, 130401 (2011). 
[12] B. Yurke, A. N. Pargellis, T. Kovacs, and D. A. Huse, 

Phys. Rev. E 47, 1525 (1993). 
[13] A. J. Bray, A. J. Briant, and D. K. Jervis, 

Phys. Rev. Lett. 84, 1503 (2000). 
[14] L. Mat hey, K. J. Giinter, J. Dalibard, and 

A. Polkovnikov, |arXiv: 11 12.1204] (2011). 
[15] J. Choi, W. J. Kwon, and Y. Shin, Phys. Rev. Lett. 108, 

035301 (2012). 
[16] J. Choi et al, New J. Phys. 14, 053013 (2012). 
[17] M. S. Heo, J. Choi, and Y. Shin, Phys. Rev. A 83, 013622 
(2011). 

[18] S. Dettmer et al, Phys. Rev. Lett. 87, 160406 (2001). 

[19] S. Manz et a/., Phys. Rev. A 81, 031610(R) (2010). 

[20] The resolution of our imaging system is about 5 /im, so 
the spectral signal at q > 1 /im" 1 is purely contributed 
from photon shot noises. We subtract the constant offset 
value at high q from the spectrum, which is typically a 
few%ofP(gi). 

[21] H. F. Talbot, Philos. Mag. 9, 401 (1836); L. Rayleigh, 
ibid. 11, 196 (1881). 

[22] J. T. Winthrop and C. R. Worthington, J. Opt. Soc. Am. 
55, 373 (1965). 

[23] A. W. Lohmann and J. A. Thomas, Appl. Opt. 29, 4337 
(1990). 

[24] M. S. Chapman et a/., Phys. Rev. A 51, R14 (1995). 
[25] A. Imambekov et a/., Phys. Rev. A 80, 033604 (2009). 

[26] In our experiment, ^/^f /6d = \/te~fjJh < 6. The focal 

depth of our imaging system is ~ 100 /im, limiting the 

expansion time t e < 40 ms. 
[27] Z. Hadzibabic, P. Kriiger, M. Cheneau, S. P. Rath, and 

J. Dalibard, New J. Phys. 10, 045006 (2008). 
[28] H. Holzmann, M. Chevallier, and W. Krauth, Europhys. 

Lett. 82, 30001 (2008). 
[29] At low temperatures, the density profile of the coherent 

part is close to a parabola shape, but the gaussian fit still 



5 



gives a good estimation for n c within a few %. 
[30] T. W. Neely, E. C. Samson, A. S. Bradley, M. J. Davis, 

and B. P. Anderson, Phys. Rev. Lett. 104, 160401 (2010). 
[31] S. J. Rooney, P. B. Blakie, B. P. Anderson, and 

A. S. Bradley, Phys. Rev. A 84, 023637 (2011). 
[32] L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattore, 



and D. M. Stamper-Kurn, Nature (London) 443, 312 
(2006). 

[33] S .Hofferberth, I. Lesanovsky, B. Fisher, T. Schumm, and 

J. Schiedmayer, Nature (London) 449, 324 (2006). 
[34] G. Roumpos et a/., PNAS 109, 6467 (2012). 



